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STEADY-STATE HEAT TRANSFER IN TRANSVERSELY HEATED POROUS MEDIA 
WITH APPLICATION TO FOCUSED SOLAR ENERGY COLLECTORS 

by Lester D. Nichols 
Lewis Research Center 

SUMMARY 

A fluid flowing in a porous medium which is being heated from a direction that is 
transverse to the fluid flow is considered. This configuration resembles the design of 
a focused solar energy collector for use in an electric power generating system. The 
effect on working fluid outlet temperature of the conductivity of the porous medium, 
both transverse and parallel to the flow is evaluated, as well as the heat transfer be- 
tween the fluid and the porous medium. 

A fluidized bed can be regarded as a porous medium with special properties. An 
analysis of a fluidized bed collectii^ concentrated solar energy to heat the working fluid 
for a heat ei^ne indicates the advantage of high thermal conductivity in the transverse 
direction and high melting temperatures of the fluidized material. 

The collector will lose heat to the environment from the absorbing surface via both 
convection and radiation. The radiation loss depends on the emissivity of the surface 
and can be controlled by use of window materials with appropriate spectral character- 
istics. The analysis indicates the effect on collector performance of the radiation loss 
when added to the always present convection loss. 


INTRODUCTION 

Solar energy can be converted into electrical energy by using it as a heat source 
for a closed thermodynamic cycle power generating system. The efficiency of such a 
system improves as the maximum temperature of the working fluid increases. This 
maximum temperature can be increased by focusing the solar energy but is limited by 
collector design. Cycles using gas or vapor as a working fluid are the least expensive. 
However, these working fluids are nearly transparent and cannot absorb the solar 
energy. In order to use these systems, the collector m\ist be designed to absorb the 



solar energy. This is commonly accomplished by an opaque surface. The heat is 
transferred from the opaque surface to the working fluid by conduction and convection. 
This step must be carefully considered, because it can introduce large temperature 
drops between the opaque surface and the working fluid, particularly at high radiant 
flux levels, and thus decrease system efficiency. Or it can introduce large heat trans- 
fer area (and consequently, large equipment size) in order to achieve adequate heat- 
transfer rates. This will tend to increase capital costs. 

A fluidized bed is an excellent candidate for heating a gas working fluid. A fluid- 
ized bed is a container of particles through which gas flows from the bottom to the top. 
The gas velocity is adjusted so the drag force on the particles is equal to the particles' 
weight, so that the particles levitate. This suspension has many of the properties of a 
fluid, hence its name. If one of the vertical walls is transparent (fig. 1) then the solar 
flux can be absorbed by the particles. Fluidized beds have very high effective conduc- 
tivity (ref. 1), which means large heat-transfer rates from the absorbing surface into 
the interior where the gas can pick up the energy. Not only that, but the effective con- 
ductivity is larger even for materials without large thermal conductivity (e.g. , ceram- 
ics) which have high melting points. Thus, the use of a fluidized bed to transfer solar 
energy to a gas may remove the limitations of either low meltii^ temperatures for 
materials with high conductivity, or low conductivity for materials with high melting 
points. 

A fluidized bed collector can be analyzed as a porous medium with through-flow. 

In the analysis the temperature of the porous medium and the temperature of the gas 
are calculated from the energy conservation equations. The three-dimensional, time 
dependent equations have been derived by Bland (ref. 2). The local velocity must also 
be specified in these equations. It is specified as proportional to the pressure gradient 
by using D'Arcy's Law. In heat-transfer cases such as this where the pressure is con- 
stant across the flow area, the velocity profile is uniform and the flow is called plug 
flow. 

The coupling between the temperature of the porous meditim and the temperature of 
the flowing gas is dependent on the coefficient of heat transfer. If the coefficient is 
sufficiently large or particles sufficiently small, then the temperatures will be equal, 
except in a small region near the gas flow entrance to accomit for a difference in the 
boundary condition (ref. 2). Needless to say, this simplifies the analysis considerably 
by eliminating one of the energy equations. 

In this particular application the heat transfer is transverse as well as parallel to 
the flow direction. This means the analysis must be two dimensional. Not only that, 
but experiments indicate that the effective conductivity of the fluidized bed is signifi- 
cantly larger in the direction parallel to the flow as compared to the direction trans- 
verse to the flow. The analysis must consider a non-isotropic effective thermal 


2 





conductivity. But only the steady state is considered. Finally, the boundary conditions 
on the porous medium for this application specify the heat flux impinging on the surface. 
They do not specify the temperature of the surface. 

Heat transfer between a flowing fluid and a porous medium has been studied inten- 
sively. Most of these studies have been one dimensional with the fluid flowing across 
the surface upon which the heat flux is applied (ref. 3). The application is transpiration 
cooling. However, there have been some two-dimensional studies. An important ap- 
plication of flow in porous media occurs in geological applications (refs. 4 and 5). The 
phenomena are complicated and the solutions are numerical. Bland (ref. 2) presented 
a derivation of the equation for a perfect gas and presented analytical solutions for some 
special cases, but none for the two-dimensional steady-state case with the heat-flux 
boundary condition. Siegel and Goldstein (ref. 6) obtained analytical solutions for two- 
dimensional cases where the fluid is flowing across the surface on which the heat flux 
is applied with either the heat flux or temperature boundary condition. They treated a 
variable thickness porous medium and used D'Arcy's law for the velocity of an incom- 
pressible gas. But they considered the case where the heat- transfer coefficient between 
the porous medium and the gas is infinite. Siegel (ref. 7) extended this work to include 
flow parallel to the heat flux svurface (i.e. , a coolant passage). However, he obtained 
solutions for only the case of a specified surface temperature and not the specified heat 
flux. Koh and Colony (ref. 8) considered the two-dimensional problem with the heat 
flux specified for plug flow in a channel of constant cross sectional area without assum- 
ing the heat- transfer coefficient between the gas and the porous media to be infinite. 
However, they did not consider nonisotropic conductivity and obtained only numerical 
solutions. 

In this report the steady-state solution to the two-dimensional problem with heat 
flux boundary conditions and a finite heat-transfer coefficient is obtained. Plug flow in 
a constant area flow passage and constant specific heat and thermal conductivity is 
assumed. The thermal conductivity can, however, be nonisotropic. Solutions are also 
obtained for limiting values of the parameters. An appropriate limiting solution is then 
used to evaluate the absorber as the heat source for a Carnot thermodynamic cycle. 
Calculations illustrate the advantage of the trade-off between absorber high melting 
point and high thermal conductivity and the effect of collector losses to the environment. 


ANALYSIS 
Governing Equations 


Consider a porous medium contained in a transparent walled container as shown in 
figure 1. This container is the radiation absorber. A flux # impinges on the trans- 
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parent wall, passes through the wall, and is absorbed by the solid material. The ab- 
sorbed energy is conducted away from the transparent wall in the porous media, and 
transferred to the gas flowing in the pores by convection. Some is lost from the porous 
medium to the environment. The heat conduction parallel to the flow is assumed to be 
different from the heat conduction perpendicular to the flow. The steady-state energy 
balance equation for the gas in plug flow is (ref. 9, p. 530, eq. 42- 1) 


(Pu)«c„ -^ = k, 




'op 


3z 


sx 


+ k, 




sz 


dx‘ 


3z" 


( 1 ) 


(The symbols are defined in appendix A. ) 

This equation is valid when the chosen coordinates coincide with the axes of prin- 
ciple thermal conductivity (ref. 10, p. 14). Such will be the case for a fluidized bed 
oriented with the z axis parallel to the body force. The conductivity of the medium in 
both directions is an effective conductivity. It depends not only on the type of medium 
but also the type of gas (ref. 11). Heat conduction by the gas is neglected since the 
effective conductivity in the medium is t 3 q>ically 10 times that of the gas (ref. 11). 
There is no net motion of the solids. The steady-state energy balance for the medium 
is (ref. 9, p. 530, eq. 42-2) 


a^t. 


sx 


+ k, 


a^t. 


Sh, 


ax^ 
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(ts - 


az" 


( 2 ) 


Boundary Conditions 


The gas temperature at the entrance (z = 0) will be given as t^. The boundary con- 
ditions specified on the porous medium is the heat flux. At the entrance (z = 0) and the 
exit (z = L), the heat is transferred from the porous medium to the gas by convection. 
The boundary condition should be (ref. 12): 
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3z 


lz=0 
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= KK - W 

dz 


^ex^^ex " ^s^ 


lz=L 


Integrating equation (1) and appljdng these boimdary conditions indicates that if 
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at z = 0 


(puLc. 


op 
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and 


(Pu)qC t 

2_E>> 1 - _e. 


‘ex 


'ex 


at z = L 


then the energy lost by the porous medium end surfaces to the gas is negligible com- 
pared with the energy picked up by the gas from the porous media volume. These in- 
equalities apply for all cases considered herein so that only a small error is introduced 
when boundary conditions are 


at, 


s 

dz 


= 0 and 


z=0 


0Z 


= 0 


z=L 


We will consider only this boundary condition in the z direction. In the x direction 
no heat will be lost from the porous medium along the insulating plane (x = 0). However, 
at the absorbing surface, (x = B) the absorbed energy will be either conducted inward or 
lost to the environment: 


at. 


sx 


0X 


+ - *env) = ^ 


lx=B 


(3) 


Nondimensional Variables and Parameters 


Introduce the following variables : 


I 






hw[t(x,z)-t^^y]- 

^w^^o ~ ^env^ " ^ 


(4) 
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" ^env^ " ^ 


and parameters 




V 


Lh, 


w 


Bh^ Lh„ 


( 5 ) 


into equations (1) and (2). The equations and boundary conditions become 


’ ~ ~ ^(^ ” ^) 

(6) 

<p(?,0) = l 

(7) 

^^{0,0 = 0 

(8) 

'HI, ?) = o 

(9) 

^^(^»o) = */^^(^>l) ~ 0 

(10) 

General Solution 



The separation of variables technique will be used to get the solutions. ^ Assume 
solutions of the form 


<P ajj(C)cos 
n=0 


( 11 ) 


“n^ 

n=0 


( 12 ) 


Equations (8) and (9) are satisfied if 


ItWs technique was used by Bland (ref. 2) to solve a time dependent problem with one space 
coordinate. His solutions are the same as herein if kg 2 aud x = UqT where t is time. 

However, he considers different boundary conditions. 
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n = 0, 1, 2, . . . 


(13) 


ua^ sin = cos 

Substitute (p and \j/ from equations (11) and (12) into equation (6): 




^bn + /3aj^)cos o-j^^ = 0 


n=0 


(14) 


V (coa^ - Iih'^ + i/Q!2bjj)cos = 0 (15) 


These equations must be valid for all hence, 

wa' - /3b„ + /3a„ = 0 
n ^ n n 

2 

cva’ + i^a„h - /j.b” = 0 
n n n n 

This set of equations possesses solutions of the form 


= A„e 

n n 




and b„ = B„e 
n n 




(16) 

(17) 


(18) 


where 


and 



^ 


(19) 


= (^ + wyji)(^^yn ■ 


( 20 ) 


Equation (20) is a cubic equation in y^^ and is the same as equation (71) in Bland 
(ref. 2). There are three roots. Let y^^^ denote the i*^ root. Then equations (11) and 
(12) can be written as 
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^ r\ 

<P = 

n=0 i=l 


( 21 ) 


and 




V y A<‘)il^e’'"’' 


cos 


( 22 ) 


n=0 i=l 


There are three boundary conditions that still must be satisfied: The two in equation (10) 
require 


i=l 

i=l 


(23) 


(24) 


and equation (7) requires 


£ 2 cos 0^1 = 1 
n=0 i=l 


(25) 


But, since the eigenfunctions are orthogonal on the interval 0 to 1, equation (25) can be 
put into the form of equations (23) and (24) by multipl 3 dng by cos and integrating 
over the interval 0 to 1 and only terms m = n remain: 
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Z 

i=l 


.(i)-. 


2 sin a 


n 


- = M 

0!„ + Sin 0!„ cos of„ 
n n n 


(26) 


Equations (23). (24), and (26) represent three equations for the three unknowns 

and Ap'. The solutions are 
n ' n 


II 


( 27 ) 


and 
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(3) 


,(1)' 

^n j 




/3_ 


(28) 


In this way the solutions for <p and \{/ are given by equations (21) and (22) with the 
Q!j^'s given in equation (13), the y^^'s given by equation (20), and the A^^^’s given in 
equations (27) and (28). Solutions are presented in appendix D for limiting values of 
some of the parameters. 

There are two numbers which are important in the application of porous medium 
collectors: the amount of heat picked up by the gas and the maximum temperatime of 
the porous medium. The heat transferred to the gas increases the average temperature 
of the gas, which is defined as 



Then the average temperature from equation (11) is 


(29) 
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(30) 


*/IV 


oo 

(t> - t- \ ' sin Q!_ 

+ *env - *0 /Lj “n 


n=0 


The maximum temperature of the porous media, = t (x = B, z = L) can be written 

from equation (12) as 


(31) 


COLLECTOR PERFORMANCE AS THE HEAT SOtJRCE FOR A CARNOT CYCLE 

The collector analyzed herein can be considered as the heat source for any closed 
thermod 3 mamic cycle. The various collector parameters and their effects on the effi- 
ciency of such a system when the closed thermodynamic cycle is a Carnot cycle (or has 
an efficiency that is a fixed fraction of a Carnot cycle) are studied. The efficiency of a 
Carnot cycle operating between the average temperature of the gas at the outlet of the 
collector and the temperature of the environment is calculated. 

The fraction of the power impinging on the collector that is collected by the gas is 


^coll 


(pu)oC B«t) - t^„^) 


$L 


(32) 


and the efficiency of a system with Carnot efficiency is 




<t) - t, 


env 


<t> 


(33) 


This can be written in terms of oj (eq. (5)) and 0 where 


0 = 


$ 


^w^env 


(34) 


as 
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t] = 


0 


*env<t> 




It is the intention here to determine the maximum efficiency possible for a collector of 
given size and of a given porous medium. Intuitively one expects this to occur at the 
maximum amount of heat that the gas can absorb. Therefore, the temperature of the 
gas and the porous medium are assumed to be the same at the exit of the collector. This 
can be achieved if L is large. Just how large will be shown later. The conductivity of 
the porous medium in the direction parallel to the gas flow is assumed to be infinite. 

This restriction will tend to lower the exit gas temperature (since the porous medium 
will be at a constant temperature in the parallel direction) but is consistent with the 
properties of a fluidized bed. 

Under these two conditions the limiting solutions are obtained in appendix D as 

* " ‘s (‘s - 


and 


( 37 ) 

where x is defined in equation (D9) and all boundary conditions except equation (3) are 
satisfied. This boundary condition relates the heat flux to the collector characteristics. 
In appendix D, C is evaluated using equation (3) when the heat loss from the collector 
to the environment is by convection only. However, at high collector temperature, heat 
loss to the environment may be dominated by radiation. In this case the boundary con- 
ditions should be 


*s I ‘ 

‘-IT "4 



(38) 


If 


K = 


env 


(39) 


the equation to determine C becomes 


Ci^x sinh X + — |(1 + C cosh x)^ - l| = 0 


(40) 



This definition of is chosen so that, when tg is only slightly different from 
equation (38) will reduce to equation (3). Equation (40) defines C which in turn defines 
tg. The average temperature of the gas leaving the collector is 


<t> 



(41) 


If all the parameters are fixed except the mass flow parameter, there is an optimum 
value of the mass flow parameter that maximizes the efficiency. This optimum occurs 
because, when the mass flow is too low, the collector transfers a very small amount of 
heat to the gas and the collector losses are a large fraction of the impinging flux. On 
the other hand, at large values of the mass flow parameter, the exit temperature of the 
gas is low and the efficiency of the Carnot cycle is very low. The calculations for the 
performance of the collector as a heat source for a thermodynamic cycle will be made 
at this optimiim value of the mass flow parameter. 


CALCULATION PROCEDURE 

It is necessary to evaluate the eigenvalues from equation (13) in order to obtain 
numerical answers from the solutions given in equations (21) and (22). These eigen- 
values were obtained using the numerical technique described in appendix B. After cal- 
culating each eigenvalue, it is necessary to obtain all three roots of the cubic equa- 
tion (20). The procedure used in the calculation, along with some limiting values, is 
shown in appendix C. The calculations were made on a desk calculator, which had been 
equipped with additional memory. 


RESULTS AND DISCUSSION 

To discuss the effect of some parameters on the average outlet gas temperature and 
the maximum porous medium temperature, let us consider certain of the parameters as 
fixed. Consider a given radiative absorber size (fixed L and B) with a specified heat- 
transfer parameter to the environment h^. We will then look at the effects of fluid- 
particle heat transfer, parallel and transverse porous medium conductivity, and the 
mass flow of gas passing through the absorber. Many studies of flow in porous media 
have assumed that the heat-transfer coefficient between the gas and the medium is 
infinite. This means that the gas and medium temperatures are equal (except in a small 
entrance region). First the condition under which this assumption is valid in so far as 
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its effect on the average gas temperature at the absorber exit is concerned is exam- 
ined. 


Mass Flow Parameter Chosen Independently 

Effects of heat- transfer coefficient . - In figure 2 the average exit gas temperature 

and maximum temperature of the porous medium are examined for three different values 

of the mass flow parameter. The transverse and parallel conductivity parameters are 

set equal. Consider first figure 2(a) where the conductivity parameters are 0. 01. 

There is a value of /3 (depending upon o)) beyond which further increases have negligible 

effect on either the maximum temperature of the porous medium or the average gas 

temperature. This value is about 8a>. Examination of figure 2(b) for conductivity 

parameters of 1. 0 reveals the same conclusion. As a matter of fact, equations (D3) and 

(D4) show the same results when p and u become infinite. And equations (Dll) and 

(D12) show the same thing when only ju is infinite. This condition can be related to a 

minimum desirable length of the absorber. The condition /3 > 8u) can be written as 

L > 8(pu)-C V/Sh for a fluidized bed of spheres, where the void fraction is talcen to 
op s 

be 0. 4. A good approximation for the gas-to-particle heat-transfer coefficient for gases 
with Prandtl number of 1, and particle Reimolds numbers greater than 1 is (ref. 13, 
p. 215) 


= — (Pu)^c_ 
s 20 ° P 

Substituting the surface to volume ratio and heat-transfer coefficient into the minimum 
length requirement 3 rields the following relation: 

L > 45 d„ 

Effect of parallel c onductivity of porous medixun . - To examine the effects of the 
parallel conductivity parameter, take, for example, the case of reasonably thick col- 
lectors in the atmosphere where ^ = 100 might be typical. A reasonably high conduc- 
tivity in the transverse direction (say copper) would 3 deM v = 100. For these param- 
eters the effect of mass flow on both the average gas temperature and the maximum 
soUd medium temperature is shown in figure 3(a) for three values of the ratio of parallel 
conductivity to wall heat loss. Notice that high average temperature requires low mass 
flow rates. Also, note that, as the ratio of parallel to transverse conductivity increases, 
the average temperature decreases (except at the extreme limits of no temperature in- 
creases or no gas flow into the absorber). Note finally that there can be a difference in 
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temperature between the maximum porous medium temperature and the average gas 
temperature (except in the limit of no mass flow into the absorber). 

Now consider an abosrber with a transverse conductivity which is much smaller 
than the previous absorber, say porous alumina. A typical value here might be v = 1. 
The same parameters are shown for this absorber in figure 3(b). For this case the 
temperature difference between the average gas temperature and the maximum porous 
medium temperature is negligible. A comparison of figures 3(a) and (b) reveals that the 
effect of the poorer solid conductivity can be compensated for by reducing the mass flow 
so that the same fluid average temperature can be maintained. The curves in figure 3 
were calculated using the general solution. However, as the parallel conductivity 
parameter ju becomes larger, the limiting form of the solution obtained in appendix D 
can be used. At ju = 100 the two solutions are indistinguishable (fig. 3). 

Effect of transverse porous medium conductivity . - In figure 4 the dimensionless 
values of the average exit gas temperature and maximum porous medium temperature 
are shown for three different values of the mass flow parameter. The heat-transfer 
parameter j8 is 100 and the transverse and parallel thermal conductivity parameters 
are equal. At low values of the transverse conductivity, a significant difference between 
the maximum and average temperature exists. However, the difference become less 
than a few percent for values of the transverse conductivity parameter greater than 
yiOct). As a matter of fact, when this condition is met, the limiting solution given in 
appendix D (eqs. (D3) and (D4)) is correct. Thus, for equal transverse and parallel 
conductivity parameters, the condition specifies a minimum value of the conductivity 
parameter for which the transverse temperature gradients can be neglected. Larger 
values are unnecessary. This implies that for a given porous material the thickness B 
need be no smaller than some minimum value (see eq. (5)). 

Effect of mass flow . - In figure 5 the dimensionless values of the average exit gas 
temperature and the maximum porous medivun temperature are shown as a function of 
mass flow parameter for several values of the transverse conductivity parameter. 

Again, the heat-transfer parameter ^ is 100. However, the parallel porous medium 
conductivity parameter pi is taken as infinite. This is a reasonable approximation for 
a fluidized bed with L/B less than 10. In this case, the limiting form of the solution 
as given in appendix D (eqs. (Dll) and (D12)) is valid and is used in making the calcula- 
tion. Once again, a minimum value of the transverse conductivity parameter for which 
the transverse temperature gradients can be neglected is found. The value is 10 and, 
in turn, specifies a minimum thickness for the porous medium when a material is 
selected. 

Up to now the mass flow parameter has been assumed to be independent. We have 
found that the average temperature of the gas leaving the absorber does not change when 
the collector length exceeds a certain value, when the collector width becomes smaller 
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than a certain value, or when the ratio of parallel to transverse conductivity exceeds a 
certain value. Some of these limiting values depend on the mass flow parameter value. 
For certain applications there may be other considerations which are important and they 
dictate that a certain mass flow parameter be used. In this case the mass flow param- 
eter Avill no longer be independent. One such application is using the collector to heat 
the working fluid for a closed thermodsmamic cycle. 


Mass Flow Parameter Dependent On Heat Flux and Transverse Conductivity Parameters 

Consider the collector as the working fluid heater for a closed thermodynamic cycle 
with Carnot efficiency or a fixed fraction of Carnot efficiency. Intuitively, one feels 
that this heater should be designed to have a minimum temperature difference between 
the porous medium and the average temperature of the gas leaving the heater. This 
impUes a collector length greater than some minimum (in the case of spherical particles 
it is 45 dg). If it is assumed that that condition is satisfied, an infinite heat-transfer 
parameter can be used to make calculations. 

The large effective conductivity of a fluidized bed makes it attractive for this appli- 
cation. The parallel conductivity in a fluidized bed is much larger than the transverse. 
For beds with length to width ratios no greater than 10, the transverse conductivity 
parameter p will be large. As long as it is larger than 10, the calculations can be 
made assuming the parameter is infinite. Thus, all calculations will be made using this 
limiting solution. 

The effect of the two remaining parameters, the transverse conductivity parameter 
V and the heat flux parameter 0, is now considered. The mass flow parameter will 
be expressed in terms of these two remaining parameters. It will be determined as be- 
ing the value that maximizes the cycle efficiency. 

The efficiency of the cycle also depends on the t 3 ?pe of loss from the collector to the 
environment. All collectors will have a convection loss which depends on their surface 
temperatures and the conditions of the environment. In addition, there may be radiation 
losses to the environment. These losses can become large as the surface temperatvire 
increases, but they can be counterracted by using a low emissivity surface. This can 
be accomplished to some degree by using a window that transmits the short wave length 
incoming flux and reflects the longer wave length outgoing flux. In this way, even though 
the collector has a high absorptivity, its emissivity may be small, and the window may 
reflect the absorber's emitted energy back. The window losses to the environment 
would then be governed by convection, even at high temperatiires. In order to consider 
the effect of this option, two cases are presented. The first is a case with convection 
only, regardless of collector temperature. The second is a case with radiation only. 
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regardless of collector temperature. However, in the second case, the emissivity of 
the radiating sur&ice is assumed to be adjusted so that at low surface temperatures, the 
radiant heat loss is numerically equal to the convective teat loss. The radiating surface 
is assumed to be at the same temperature as the collector surface. 

The efficiency is calculated for a Carnot cycle, with losses to the environment by 
convection only, as a function of mass flow parameter (fig. 6). There are two different 
values of transverse conductivity parameter and three different values of the heat flux 
parameter. Note that there is a value of mass flow parameter that maximizes the cycle 
efficiency and that its value depends on the remaining two parameters. All subsequent 
calculations of collector performance assume that the mass flow parameter is optimum. 
The effects of varying the two independent parameters on the cycle efficiency and the 
optimum mass-flow-rate parameter are investigated. And these quantities are plotted 
as functions of the ratio of the maximum temperature of the porous medium to the en- 
vironment temperature. This method of presentation is then useful when comparing the 
two different types of collector losses for the same porous medium. 

The calculations for the efficiency are shown in figure 7 for the two heat loss con- 
ditions and illustrate the expected results that as the ratio of medium to environment 
temperature increases and the conduction parameter increases the efficiency increases. 
The efficiency approaches a maximum value monotonically for the convection condition 
as the flux parameter increases for a fixed value of the conductivity parameter. How- 
ever, for the radiation- loss case (fig. 7(b)), there is a maximum efficiency as the flux 
parameters increase for fixed values of the conductivity parameter. 

As one expects intuitively, at a given value of the heat flux parameter, the effi- 
ciency for the convection boundary condition is better than for the radiation condition - 
regardless of the value of the conductivity parameter. However, for a given ratio of 
medium to environment temperature, the efficiency with the radiation- loss condition is 
higher than with the convection- loss condition, if the conduction parameter is large 
enough. These higher efficiencies, however, require much larger heat flux parameters. 
Comparing figures 8(a) and (b) shows that much larger mass flow parameters are also 
required for the radiation loss. Thus, even though the heat losses per unit receiver 
area are larger (at a fixed porous medium temperature) with the radiation- loss condi- 
tion, the total area required is smaller. This is true because of the higher mass flow 
rates required at the high heat flux values. The result is that a larger fraction of the 
incoming flux is transferred to the gas. 

To achieve these high efficiencies, a very large value of the heat conduction param- 
eter is also required. Because of its high effective conductivity, a fluidized bed of 
particles can attain the needed value in a bed which is much thicker than a packed bed 
using copper. Also, if the particles are made of a material with high melting tempera- 
ture, the efficiency of the system can also be improved by operating at a higher tem- 
perature. 
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PERFORMANCE OF A SOLAR COLLECTOR 


The performance characteristics of solar collectors using porous copper and fluid- 
ized alumina are compared in figure 9. The length L to thickness ratio B for all 
collectors is 10. The two heat- loss boundary conditions are shown for one temperature 
ratio (3) for copper and the radiation- loss condition for three temperature ratios (3, 4, 

5) for the fluidized alumina. Results are shown for three thicknesses: 0. 1, 1, and 
10 meters. The conductivity is 70 watts per meter per kelvin for sintered porous cop- 
per (ref. 14) and 800 watts per meter per kelvin for the fluidized alumina (ref. 1). The 
heat-transfer coefficient is 9. 3 watts per square meter per kelvin for convection 
(ref. 15, p. 529), and 8. 1 watts per square meter per kelvin for radiation with an en- 
vironment temperature of 330 K. The first is the natural convection coefficient for a 
long vertical flat plate in air, and the second is for a black absorber at the environment 
temperature. 

Some trends and conclusions from these comparisons are as follows. Figure 9(a) 
shows that absorber thicknesses on the order of a meter for the fluidized alumina pro- 
vides the same efficiency as copper with thicknesses on the order of 0. 1 meter. This 
implies larger mass flow and lower capital cost. At a temperatiire ratio of 3 and a 
thickness of 0. 1 meter for copper, the efficiency for both the convection and radiation 
loss cases are about the same; however, at larger thickness the efficiency of the system 
with the radiation boundary condition drops off rapidly. Generally, a system has higher 
efficiency if the temperature ratio is high and the thickness small. 

The heat flux which must be applied to achieve the efficiencies is shown in fig- 
ure 9(b). For a radiant flux of 640 kilowatts per square meter (used for design of the 
system in ref. 16) the fluidized bed could operate satisfactorily with a temperature ratio 
of about 3. 5 and a thickness of about 1 meter. A Carnot cycle would have an efficiency 
of about 55 percent. If an actual efficiency of 0. 7 of Carnot could be achieved, then the 
efficiency could be 38 percent. These results are comparable with those of reference 17 
where it is suggested that a liquid metal might be required for such large heat fluxes. 

The mass flow rate parameter is shown in figure 9(c). Immediately, one can see 
that approximately 10 times as much mass flow is required for the radiation- loss con- 
dition for the convection- loss condition for the porous copper. Also, for the heat flux 
of 640 kilowatts per square meter one can calculate that for air at a few atmospheres 
pressure a fluid velocity of a few meters per second is required for a fluidized bed. 

This is within the acceptable limits for bed design. Thus a fluidized bed may circumvent 
the need for a liquid- metal working fluid at high heat flux conditions. 
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CONCLUSIONS 


The high effective conductivity of a fluidized bed makes it an attractive candidate as 
an absorber of thermal radiation for the purpose of heating an otherwise transparent 
fluid. The particle to gas heat-transfer coefficient is sufficiently large so that the tem- 
perature difference between the gas and the solid can be negligible. The longitudinal 
conductivity is much larger than the transverse conductivity, but the effect on the gas 
temperature can be coimterracted by adjusting the other parameters. Finally, prelim- 
inary calculations for electric power systems on Earth indicate that the fluidized bed 
absorber could provide the required performance. Of course, some factors not con- 
sidered in this ideal calculation could have a detrimental effect and would have to be 
considered for any actual applications. These include, among others, the effect of 
pressure drop in the collector on systems performance and the deterioration of the 
optical transmission through the transparent wall because of particle scratches or other 
contamination of the wall. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 18, 1976, 

506-24. 
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APPENDIX A 


SYMBOI^ 

A expansion coefficient defined in eq. (18) 

a expansion coefficient defined in eq. (11) 

B expansion coefficient defined in eq. (18) 

b expansion coefficient defined in eq. (12) 

C constant of integration defined in eq. (56) 

c specific heat 

D determinant of coefficients defined in eq. (28) 
d diameter 

h heat-transfer coefficient 

k thermal conductivity 

L length of porous media bed 

M expansion coefficient defined in eq. (26) 

S heat-transfer contact area between porous medium and gas 
t temperature 

u gas velocity 

V volume of porous media bed 
X, z coordinates shown in fig. 1 

a Eigenvalues defined in eq, (13) 

/3 heat-transfer coefficient parameter defined in eq. (15) 

€ emissivity 

? fraction of distance from entrance to exit of porous media 

Tj Carnot cycle efficiency 

0 heat flux parameter defined in eq. (53) 

M parallel thermal conductivity parameter defined in eq. (5) 

V transverse thermal conductivity parameter defined in eq. (5) 

4 fraction of distance from insulating wall to absorbing 

p gas density 


a Stefan- Boltzmann constant 

$ heat flux absorbed by porous media bed 
<p dimensionless temperature of gas defined in eq. (4) 

X parameter defined in eq. (4) 

xj/ dimensionless temperature of porous media defined in eq. (4) 

(t) mass-flow parameter defined in eq. (5) 

Subscripts; 

coll collector 

env environment 

ex exit 

g gas 

max maximiun 

n, m summation indices 

o entrance 

s solid porous media 

sx transverse 

sz parallel 

w wall 

Superscripts: 

(i) summation index 

~ parameters and variables defined in appendix D for solution with no heat lost 

from absorbing surfeice 

*, first and second derivatives with respect to ? 
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APPENDIX B 


EVALUATION OF EIGENVALUES 
It is necessary to find from equation (13): 

"'“n* = -;7 '=“* ”n 
“n 


The Newton-Raphson technique can be used if a good guess for the initial trial can be 
made. One which works can be found by replacing 


cot with 



This form has the correct zeroes and infinities. Then, the positive root is 






umT + 2vn + {um - 1 ) 
2v 


(Bl) 


The convergence of the Newton-Raphson technique is quicker if the logarithm of v is 
taken: 


d]n u 1 

Q!„ Sin a„ cos 
n n n n 

Approximating the derivative by a difference and evaluating it at result in 


JL- 1 

t . t t 

sin a cos 
n n n 


In - In 


cot at 


a 


n 


n n 


(B2) 


(B3) 
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Thus, 


a 


n 



cot a^, 
In 1 / - In 



t . t t 

q;„ sin Q!„ cos 0 !„ 

n n n 


(B4) 


After each trial is taken as the oi^. Convergence to any desired accuracy can be 
achieved. Obvious limiting forms are 


Um a 

I /-0 


n 



lim 


a 


n 


HTT 


n :?t 0 


0 J 


(B5) 
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APPENDIX C 


ROOTS OF CUBIC EQUATION 


Equation (20) can be rewritten as 



//3 + 


n 




Zj) 3u,\9^2 m 



where 




A convenient form to write the roots is achieved (ref. 18) by letting 


cos <p 


g 4 ,3 - 

60.(9 ^2 II j 

/ 2\^^^ 

1 


Then 



(Cl) 


(C2) 


(C3) 


(C4) 
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APPENDIX D 


UMITING SOLUTIONS 

Infinite Conductivity of Porous Medium (f— ; M-*°° ) 

Equations (1) and (2) can be combined to yield an equation for the temperature of 
the gas: 


(pu)oCp 


3z 



(ts - t) 


The energy balance on the porous medium requires 

B(pu)oCp[t(D - t(o)] = [$ - h^(tg - tg^^)]L 


(Dl) 


(D2) 


where t(o) = t^. 

The solutions for the temperatures of the gas and the solid, using the parameters 
defined in equations (4) and (5) are 


and 




(D3) 


^ = 


1 - e~^^^ 
1 + 1 - 


(D4) 


Infinite Conductivity of Porous Medium in the Direction Parallel to Flow (m— “ ) 

In this case the temperature of the porous medium is uniform in the z direction. 
Equations (1) and (2) can again be combined to provide a local heat balance for the fluid. 
(See eq. (Dl). ) However, in the porous medium the heat flux in the z direction is 
determined by the boimdary condition. This condition can be introduced by integrating 
equation (1) over z. 
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0^t 


^dz.k 


sz 


at 


z=L 


s 

3z 


(D5) 


iz=0/ 


Again, we require that no heat be conducted from either end of the porous medium, even 
as becomes infinite. Also, as k ^ becomes infinite t^ becomes a constant in 
the z direction. If the order of integration and differentiation is exchanged, then equa- 
tion (D5) becomes 


(pu)QCp[t(x, L) - t(x, 0)] = kgjjL (D6) 

dx^ 

Since t_ does not depend on z, equation (Dl) can be integrated: 

t = tg + (ts - t^)e-^^/‘^ (D7) 

This equation can be used to substitute for t in equation (D6), which can then be solved. 
The solution with dt_/dx = 0 at x = 0 is 

O 

ts =yi + Ccoshx4) (D8) 


where 


x2 H ii! (1 - 

V 


(D9) 


The constant C must be determined from the boundary condition oii the flux absorbing 
surface, as given in equation (3). In this case C is 


C = 


* V^env - *q) 
t^h^(cosh x+ v\ sinh x) 


(DIO) 


The solution can be written in terms of (p and xp using equation (4) as 

= 1 - Xl(l - e 

v\ sinh X + cosh x 
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^( 4 ) = 1 

v\ sinh X + cosh x 


(D12) 


Infinite Conductivity of Porous Medium in Direction Transverse to Flow (t>— «>) 

In this case the temperature of the porous medium is \iniform in the x direction. 
Equations (1) and (2) can again be combined to provide a local heat balance for the fluid: 


— (t, - t) 


Sh, 


'op 


dz 


(Dl) 


The heat balance for the porous media must be global in the x 


direction: 


(pu)oCp 




(D13) 


As k_^ becomes infinite the integrals can be evaluated. Substituting the variables de- 
sx 

fined in equations (4) and (15) into equation (41) results in 


The solution to these equations are of the form 

3 

<P(?) = ^^ A 
i=l 



(D14) 


(D15) 
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where the 



are given as the roots to 



and 


[ ( 3 ) 

E (2) 

/o 


■N 




with 






(D16) 


(D17) 







L^o 


( 1 ) . y(S)' 







wy, 


(3)irv(2) . 


O J 






(D18) 


2 

This is equivalent to setting va^ = 1 and va^ = 0, where n = 1, 2, . . . in equation (13). 
In this case = 1 and = 0, where n = 1,2, . . . 
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Zero Heat Conductivity of Porous Media in Direction Parallel to Flow (m— 0) 


In this case /n = 0 and the adiabatic boundary conditions at the entrance and exit to 
the porous media are satisfied. Then 





and 


(D19) 


OO 



(Bva^O 

M exp 2 cos a i 


(D20) 


..y 


M 


n 


exp - 


^ + va 


n 


_o;(/3 + i/aj)] 


cos 


(D21) 


Infinite Heat-Transfer Coefficient between Porous Media and Fluid (;3— «>) 


In this case /3 — <» . The mathematical difficulty occurs because the gas and porous 
medium temperature must be equal but the gas temperature of the entrance is arbitrary. 
This difficulty is handled by allowing the derivative to be infinite at one point - the en- 
trance - so that the inlet gas temperature can be still used as a boundary condition 
(ref. 2). The for this case are given in appendix C for ^ 0; 


<P =rf/ = 


^n< 




- y,! e 


i3)Jn 


(3) 


.,( 1 ) 


-Vcos 


(D22) 


n=0 


•y'~ 'e 
^n 


(iL^n 
- ^n ® 
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No Heat Loss from Heated Surface (h^ —0) 

In this case h = 0. The parameters used so far are defined with h in the de- 

vv w 

nominator. New parameters must be defined: 


CO 

CO = — 
V 



(D21) 


The dependent variables used so far are defined so that they approach 1 as h becomes 
zero. New dependent variables must be defined: 


cp = 


(t - 


. _ <‘s - *q)^sx 


The new variables can be expressed in terms of the old: 


9 = 


P = 



^w^^o ~ ^env^ 
$ 

^w^^o ~ ^env^ 
$ 



The solutions when h^ approaches zero are 

^ 4 (1 - 

I/— 0 V 


5^ = Ijm 4 (1 - »/^) 

v—0 V 


(D22) 


(D23) 


(D24) 


(D25) 

(D26) 


while holding u), p, and /3 constant. To take these limits, substitute equation (21), 
(22), and (25) into equations (D25) and (D26): 


30 




4/ = lim < 


" 3 
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.(i) 


u 




r<‘>CI 


COS Q!q| 


+ 




^ + ojy ' 


(i) 

n 


^(iV 

/n ^ 


cos 


(D28) 


The term n = 0 has been taken out of the summation because it must be handled sepa- 
rately in the Umitii^ process. 

Consider first the terms for n = 0. From equation (13) one has 


a = 0 
o 

and from equation (20) 



u—0 V w 


From equations (26) and (28) 


(D29) 


M„ = l 
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fJ2) (1)1 


Jo Jo 
e - e 


Substitutii^ these values into equations (27) makes it possible to take the limits 


- . 

lim — _k 
iT-O V Z 


-J2) 


1 - e 




\(2) AD 


lim 

I/— 0 V 


3- 

o __ /3 


.(1) 


rs^ 

Cx) 


.(2) 


^o L' 




/3 + ojy, 


?>] 


,(2) Jl) 

le’^o 


,( 3 ) 


lim 


lim 
v^O V 


v— 0 V 

A?) 


1 - e 




r«i«i 

a; 


1 r; e 


rv/ 

/3 0) 


(D30) 


With these values and making use of the equations for y^^\ y^^\ equations (D27) and 
(D28) may be written: 
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n=l i=l 


The terms for n greater than zero offer no special difficulty if one evaluates from 
equation (26) and from equation (13) as v becomes zero: 

lim lim a = njr 

V— 0 V 2„2 

n TT 

With this evaluation of the and the defined as the three roots of 

M3Vn = (^+2y„)(Mr^nV) (D34) 

and the defined in equations (27) and (38), then the limits in equations (D31) and 
(D32) can be taken. The solutions given in those equations are the anal 3 rtical solutions 
to the problem solved numerically by Koh and Colony (ref. 8). The solution in 
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equation (D28) converges faster if one makes use of the fact that 



since from equation (26): 



u U--0 V 2 2 

n T! 


i-1 
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Gas superficial 
mass flow, (pui^ 


Figure 1. - Schematic diagram of porous media absorber. 


Maximum porous medium temperature, Umax” V/0*^ V ^^env‘V]* <J*niension!ess 


Gas flow rate 
parameter, 



Porous medium togas heat-transfer coefficient, p - BShp/Vh^ 
(b) Transverse and parallel conductivity parameters, 1. 0. 


Figure 2. - Effect of particle to gas heat -transfer coefficient parameter on maximum temperature of porous media 
and average gas temperature. 





Maximum porous medium temperature, (t,^ 3 ^ '^•'env ’*o0- dimensionless 





Average exit gas temperature, Average exit gas temperature, 

«t> -to)/[($/hw) +<tgpy-tQ)], dimensionless «t>- tQ)/[($/h^) +ftgpy-tQ)], dimensionless 

anct and 

Maximum porous medium temperature, Maximum porous medium temperature, 

max *<‘env'V]- ‘•'mensionless ftmax‘V/C®V ^#env'V]- d'mensionless 



Figure 4. - Effect of transverse and parallel conductivity parameter on maximum temperature of porous media and 
average gas temperature. Particle to gas heat-transfer coefficient parameter, 100; transverse and parallel 
thermal conductivity parameters are equal. 


^^max ■ *oVD^^ * *^env " ^oG 


«t> - tQ)/[f^/h^l + (tgpy - tp)] Transverse thermal 

mi +u)) conductivity 



Figure 5. - Effect of gas flow rate parameter on maximum temperature of porous media and average exit gas temper- 
ature. Parallel conductivity parameter Is infinite; porous media to gas heat transfer coefficient parameter, 100. 
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Carnot cycle efficiency, 





Transverse thermal 
conductivity 



(a) All collector heat loss to environment by convection. 



Maximum surface temperature to environment temperature ratio, t^^ax^^env 


(b) All collector heat loss to environment by radiation only (except for curve marked "convection, " which 
corresponds to collector losses by convection and with infinite transverse conductivity. The 0 and 
V values do not apply to this curve. ) 

Rgure 7. - Carnot cycle efficiency as function of maximum surface temperature to environment temperature ratio. 




Optimum gas flow rate parameter. 


(a) Collector losses to environment by convection only. 



Maximum surface temperature to environment temperature ratio, t^^ax^^env 

(b) Collector losses to environment by radiation only. 

Figure 8. - Effect of surface temperature to environment temperature ratio on optimum mass flow 
parameter. 



Carnot cycle efficiency, 


Type of bed 


Porous copper 
Porous copper 
I I Fluidized alumina 


Boundary 

condition 

Convection 

Radiation 

Convection 


Numbers on top of bars are thickness, B. 
in meters 



Ratio of surface to environmental temperature, 

■ Comparison of performance parameters for collectors with length to thickness ratios of 10. 
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